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Abstract 

The DIS process on nuclei is considered in the framework of LO BK equation with local impact parameter 
dependence. The initial conditions of GBW type, found in [3], was used in the solution of BK equation. 
Integrated gluon density function and F2 nucleus structure function for different nuclei are calculated. Ob- 
tained results are compared with the different parameterizations of integrated gluon density function from 
[SI [SI El IB]- The anomalous dimensions and saturation scales for different nuclei are calculated at different 
energies. The expressions of the functional form of saturation scale are obtained for the proton and different 
nuclei as functions of the impact parameter and energy. 



1 Introduction 

The phenomenological applications of the QCD BFKL Pomeron, [JJ, and of the system of interacting Pomerons, 
[9j [10], could be investigated in the framework of the BK equation j2] [3|. This framework was used in the 
calculations of the amplitudes of the different scattering processes, see [TTJ [TJl [T31 [T3J [JJ3 [TB] for example and 
references therein. The same principals of the saturation physics were applied in the different phenomenological 
and CGC type models as well, see[I3[ia[Tj3[ia[IQ[21[23[Ml^ In these approaches, 

inspired by BK equation, the impact parameter dependence of the solution of the equation was treated only 
approximately, see for example [53[[3S], whereas the phenomenological models, such as GBW model [T71 [TS] . 
treating the impact parameter dependence neglect the evolution of the amplitude with rapidity. There are also 
the studies of the impact parameter dependence of gluon structure function together with DGLAP evolution 
of the function considered in the papers [28] [19l [29] [32] and calculations of the unintegrated gluon density 
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function in the framework of modified BK equation with the factorized dependence of impact parameter and 
momentum in initial conditions in 33, 34, 35] . Nevertheless, the investigation of the BK equation with the new 
type of the initial conditions with the non-factorized impact parameter and momentum dependence is a task 
which requires a deep attention due the importance of the BK equation in the different scattering processes 
with nuclei involved. 

In the paper [3] a first step towards the including of the impact parameter dependence into the rapidity 
evolution of the scattering amplitude through initial conditions was made. It must be stressed, that introducing 
local impact parameter dependence into evolution equation trough initial condition we miss the precise treatment 
of transverse position in the evolution kernel, that does not allow to calculate the contribution of the Pomeron 
loops into the amplitude, for example. Nevertheless, considered semi-classical framework with local impact 
parameter dependence is justified due the following observation. For the scattering on the target, which size 
is large, the momentum transfer of the Pomeron line originating from the external particle is bounded by the 
form-factors of the sources, and, therefore, it is small for the nuclei (and possibly proton) targets. Therefore, 
when we consider only "net" diagram structure, i.e. the semi-classical approximation to the problem, and when 
we do not account the Pomeron loops contribution, the constraint on the momentum transfer of the Pomeron 
line is imposed on the all Pomeron lines that justify our approximation, see more details in |12) and references 
therein. 

The DIS process on the proton, which was considered in the paper j4], allowed to find a initial condition 
for the BK evolution by the fitting of the F^ function data. In the present studies we use the same as in [3] 
calculation procedure solving LO BK equation in each point of the impact parameter space, using methods 
developed in 12 . The difference with the calculations from [3] is that now we consider different nuclei as a 
target instead the proton target in the paper [4], and, therefore, the impact parameter profile of the proton we 
change on the Wood-Saxon parameterization of the nuclei density profile. All other parameters of the initial 
profile for BK equation we take the same, in spite of the differences of the processes. The applicability of this 
assumption we will discuss latter. As it will be shown, the change of only impact parameter profile in initial 
conditions will lead to the results for integrated gluon density function similar to the results of the different 
calculations of the same function in [SJ [Tj |5] , justifying this minimal modifications of the initial conditions 
for DIS on nuclei. 

The comparison of our results with the results of other calculations is based on the fact, that due the lack 
of the high energy DIS process data for the nuclei targets, we could not perform the fitting of the data in order 
to determine all parameters for initial conditions as it was done in g]. More of that, we also could not use the 
low energy data for this purpose, because the small Bjorken x evolution begins at the energies larger than the 
energies for which the nuclei data are known. Fortunately, the existing models of the integrated gluon density, 
[5J [7J [5] , allow to extrapolate the integrated gluon density function till very small values of x, x « 10~ 7 . It 
gives to us a possibility to check how good (or bad) our calculations are and, therefore, how good (or bad) our 
initial conditions are. To our surprise, curves for the integrated gluon density of different nuclei obtaining in our 
framework are in the "window" of the curves from the [SJ[SJ[71[S] which were calculated in different approaches. 
This result, as we underlined previously, justify the use of "minimal changed" initial conditions for the nuclei 
in our calculations. 

The knowledge of the unintegrated gluon density function, i.e. solution of BK equation, allows to find 
other functions and parameters connected with the DIS process. We calculate the integrated gluon density 
function and F% structure function of the DIS process on nuclei, or, more precisely the proton-nuclei ratio of 
these functions. In this case we find a parameterization of the nuclei integrated gluon density and F% functions 
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in the terms of the corresponding proton functions and number of nucleons in the nuclei. We also calculate a 
anomalous dimension of the integrated gluon density function, similarly to the calculations of [3B] , and find a 
saturation momenta of the process as a function of impact parameter space and energy. 

This paper is organized as follows. In the next section we describe a formalism of the calculations. In Sec. 3 
we present obtained results for the integrated gluon density function. In Sec. 4 we consider the calculations of 
the F2 structure function and in Sec. 5 we present results for the anomalous dimensions of the integrated gluon 
density function of the DIS process on the proton and on the nuclei. The saturation scale calculations for the 
problem of DIS process on the proton and DIS process on the nuclei are presented in the Section 6. Section 7 
is a conclusion of the paper. 



2 The low-x structure function in the momentum representation 



In this section we shortly remind the main formulae used in the following calculations (see also [4]). The 
unintegrated gluon density function f(x, k 2 ,b) of the DIS process we find solving BK equation for the each 
point in the impact parameter space: 



a 2 



f(a 2 ,b)-f(k 2 ,b) , f(k 2 ,b) 



[4a 4 + /fc 4 ]* 

L % ><•*• " I * + ><*» L % log (f) "«••» 



where we introduced the rapidity variable y = log(l/x). In Eq[T]we assumed, that the evolution is local in the 
transverse plane, i.e. impact parameter dependence of f(y, k 2 , b) appear only throw the initial condition for the 
f(y,k 2 ,b) function 

f(y = y ,k 2 ,b) = f in (k 2 ,b) (2) 

In order to exclude part of ambiguities in the solution of BK equation which arise due the absence of the NLO 
corrections, we perform the following substitute in the equation 



obtaining 



f(y,k 2 ,b) 



N r . 



f(y,k 2 ,b)a s (k 2 ) f(y,k 2 ,b) 



(3) 



dyf(y,k 2 ,b) = ^k 



a 2 



f(a 2 ,b)-f(k 2 ,b) , f(k 2 ,b) 



k 2 \ 



[4a 4 + fc 4 ]* 
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(4) 



where y — a s y. The value of a s is a constant in the LO approximation and we consider a s as parameter of the 
model, which we borrow from the fit of DIS data performed in [1] . 

The impact factors, used latter in calculations of F% structure function, are usual impact factors of the 
problem with the three light quarks flavors of equal mass included. They have the following form (see [5S] for 
example) : 



<f> L (k,m 2 q ) = 32na 



3 r i 
2 



q=l 



dpdrj 



k 2 V (l-r 1 )p 2 (l-p) 



2 r,2 



( Q 2 p(l - p) + k 2 T](l -T)) + m 2 ) ( Q 2 p{l - p) + ml) 



(5) 
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and 

dpdi] ■ 



9=1 - 70 



fc 2 Q 2 (p 2 + (l- P ) 2 ) P (l - P ) (T7 2 + (1 -r?) 2 ) + fc 2 (^ + (1-^2) TO 2 +4p (l_ /9 ) ?? ( 1 _ ?? ) m 2 

( Q 2 p(l -p) + k 2 r)(l -T]) + to 2 ) ( Q 2 p(l - p) + to 2 ) U 

Due the including light quark masses in the calculations, the rapidity variable y (Bjorken x) in BK equation is 
also modified, see details in [T71 [TB]. For each fixed rapidity y of the process, the value rapidity taken in BK 
equation is changed 

4 TO 2 

y^y- ln(l + (7) 

The form of the function f{y,k 2 ,b) at initial rapidity, i.e. initial condition for the BK equation Eq|H has 
been borrowed from the form of GBW ansatz, [171 118) . with introduced impact parameter dependence. In the 
paper [4] the following form of the initial conditions for the DIS process on the proton was used 

f P roton(y = yo, k 2 M = -^ Q|py«*P(- ^| ) » (8) 
where a saturation scale as a function of the impact parameter is defined as following 

The S(b) function here is the proton impact parameter profile with the usual normalization properties 

S(b)d 2 b = 1 (10) 

and C is a numerical coefficient which defines a value of the saturation scale at zero impact parameter and 
initial rapidity through the proton radius 

Ql( b = 0) = ^. (11) 
The generalization of this initial condition for the case of nuclei is straightforward 

fnucleusiV = yo,k*,b) = ^^A_ e x P (-^), (12) 

where 

QUb) - ^ (13) 
with the Wood-Saxon nuclei profile function S(b) for the nucleus with A nucleons which is defined as usual 

S (b) = — 1 f°° & fl41 

4tt R\ + a 2 7T 2 R A 1 + cxp / ^Ra+VE+H^j 

with the parameters 

R A = 5.7 A 1/3 GeV- 1 , a = 2.725 GeV 1 . (15) 

It must be underlined, that the introduced Wood-Saxon nuclei profile Eq. (IT41 is a attempt of the generalized 
description of the all nuclei density profiles in one formula. In reality the two parameter Fermi model expression 
for nuclei density Eq. (fT4)l is mostly applied for heavy nuclei and it is not so correct for the light ones, see [40] . 
Therefore, the precise consideration of the processes with the light nuclei involved need a introduction of the 
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yo (xo) 


rip [L^ev ) 


c 




fflq \\srev J 


3.1 (0.045) 


7.9 


0.0855 


0.108 


0.008 



Table 1: The parameters of the model. 



different from the expression Eq. (1141) nuclei density profiles and we are not consider this particular task in the 
paper. 

The C parameter in Eq. (|13[) . as well as the value of a s , value of initial rapidity ya (xo) and masses of 
three light quarks we take the same for the both cases of DIS processes on the proton and nuclei, see TabkQ] 
There is a question, is this correct approach to use the same values of the parameters for both processes with 
the different targets. In principle, due the larger parton densities in the DIS processes on the nuclei, we could 
expect the differences in the values of a s and values of yo for these two DIS processes. Nevertheless, we must 
remember that we perform the calculation in the LO and running coupling effects are not directly included in 
our calculation scheme. We also will see, that a difference in the saturation momenta of the proton and different 
nuclei is not so large, and , therefore, this difference could not set up the large differences in value of a s for both 
processes in the LO calculations. The justification of the use of the same values of a s and the correctness of used 
parameters will be shown in the next section, where surprising coincide of our calculations of the integrated 
gluon density function with the existing parameterizations of the same function will be demonstrated. The 
same values of the C parameter and value of initial rapidity, from which the evolution equation is valid, may 
be explained by the same observations. As it seems, these values are more or less universal in the leading order 
of calculations and shows the physical boundary for the small x evolution independently from the considered 
target in DIS process. 

Last remark which concerns the presented calculations is about the kinematic range of the considered 
parameters. How it was shown in [4], in the present framework the HERA data at Q 2 < 1 GeV 2 could not 
be described. Therefore, in our calculations we limit the kinematic range of the processes by the x < xq and 
Q 2 > Qq constraints. Unfortunately it means, that the main bulk of the data of the DIS on nuclei, [4"T1 142) . 
is outside of the range aof the applicability of the model, see also remarks above about the form of the nuclei 
density profile. Therefore, in the next section, we perform the comparison of our results only with the existing 
parameterizations of the integrated gluon density function , without an introducing of the experimental data 
fit as it was done in [3]. 



3 Integrated gluon density function 

In j3] was mentioned, that the definition of the integrated gluon density function xG(x,Q 2 ) through the 
f(y,k 2 ,b) function of Eqf3]has ambiguities related with the coupling constant a s . Indeed, by definition 

xG(x,Q 2 ) = J ' d 2 b J Q ^!/(x,fc 2 ,6) (16) 
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or, in terms of f(y, fc 2 , b) function 

xG(x,Q) = jdbj ^r^j^y- (17) 

We see from Eq|T7|that correct determination of the xG(x, Q 2 ) function in terms of the f(x, k 2 , b) function must 
include the integration over the running coupling, whereas our calculations scheme includes only LO corrections. 
In LO approximation the Ea ll7l may be redefined as 

xG{x,Q 2 ) = f d 2 b f Q ^-f(x,k 2 ,b) (18) 

Qs(< Qproton >) J J k 

where the value a s (< Qp ro t on >) could be considered as a some parameter, which is not necessary the same as 
in TablefT] Similarly, the integrated gluon density function for a nucleus we could write in the following form 

xG N (x,Q 2 ) = - 1 — f d 2 b f Q ^-f(x,k 2 ,b), (19) 

a s\< W nucleus > ) J J K 

with the f(y,k 2 ,b) as a solution of corresponding BK equation for DIS process on the nucleus. Our main 
assumption, which we could check only post priori, is that in LO we could take 

"s(< Qlroton >) ~ a s{< Qnucleus >) ( 20 ) 

and, therefore, that the ratio of the integrated gluon densities of the nucleus with A nucleons and the proton 

xG N (x,Q 2 ) 
AxG{x,Q 2 ) 

does not depend on values of a s . The first check of the approximation made is simple. We calculate R A for the 
following nuclei: gold (A=197), neodymium (A=150), zinc (A=70) and neon (A=20) at 

Q 2 = 2.5, 12, 60, 120 GeV 2 . 

Obtained results, see FigfTJ we compare with the calculations of the ratio from the [6J [7J [8] , see Figj2j Figj3] 
As it seems from the FigJ2]-Figl21 in general our results for R A are in the range defined by curves obtained by 
the parameterization [HI E] for R A ratio, and clearely more differ from the [5, 8 parameterizations of the same 
ratio. Surprisingly, obtained in absolutely different framework our results somehow interpolate between [7] 
parameterization of the ratio and stay in the "window" defined at low x by "extremal" parameterizations [5j[8]. 
The closeness of all curves at the initial point of small x evolution, (x = 0.045), shows that we indeed matched 
the small x evolution of BK equation with the DGLAP framework of [5J |H1 [Tj |S] in this point. Therefore, this 
coincidence between the curves from different calculation frameworks indeed justifies the form of used initial 
conditions for the BK equation and the assumption Eq[20] It is also clear, that all parameterizations [6j [3 18] 
are based on the low energy data, whereas the high energy parts of the curves are the extrapolation of the 
established formulae in the region of small x. This is a explanation for the large differences between the curves 
from different parameterization in the region of small x in the FigJ5J 

Another question is about the parameterization of the nuclei gluon density xGn(x, Q 2 ) function in the form 

xG N (x,Q 2 ) = A aA xG(x,Q 2 ) (22) 

where is the coefficient a a (do not be confused with the coupling constant a s ) determines the "power" of the 
shadowing for each nucleus at different Q 2 . In spite of the parameterization based on the data fitting, in our 
approach this coefficient is calculable, see also [23] for the similar calculations. The results of the calculation of 
the a a coefficient is presented in the FigJU From Fig 0] it is clear, that the obtained shadowing is weaker than 
usually obtained in the framework of BK equation, see again paper [35] for example. 
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Figure 1: The ratio R A (x,Q 2 ) for different nuclei at different Q 2 . 
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Figure 2: The ratio R A (x,Q 2 ) of integrated gluon density functions: the present model results (solid lines), the 
results from £5^ (dashed lines), the results from [6] (dotted lines lines) and the results from fjj (dashed- dotted 
lines). The numbers in brackets show the added number to the real value of R A . 
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Figure 3: The ratio R A (x,Q 2 ) of integrated gluon density functions for the DIS on the gold at Q 2 = 
60,120 GeV~ 2 : the present model results (solid lines), the results from [5] (dashed lines), the results from 
[6] (dotted lines lines), the results from J^jj (dashed- dotted lines) and the results from |§J/ (stars). 

4 F2 structure function 

Using usual definition of F 2 structure function 

F ^ Q2) = ^J dH Jj^ ( < ] + ) (23) 

we obtain the same expression in terms of f(y, k 2 , b) unintegrated gluon density function 

P^Q 2 ) = S^J dH J ^^^(M*,mS) + *L(*,m2)). (24) 

The expressions for the impact factors in Ea l231 -Ea l24l are given in EqJSjEqJB] There the a s coupling constant 
is excluded from the expression comparing to the usual definition of the impact factors, that allows to write 
Eq l24l in the way where formally a s is not appeared in the expression. As in the previous case of integrated 
gluon density function, the main object of our interest is a ratio of the nucleus and proton structure functions 



pi _ F ™{x, Q 2 ) 

Rf2 ~ AF 2 (x, Q2V (25) 



The results for this ratio are presented in the FigJS] Comparing the result of the FigQ]with the results presented 
in the FigJS] it is easy to see, that in contrary to the xGn(x, Q 2 ) function the structure function F 2 n(x, Q 2 ) 
does not proportional to the AF 2 (x, Q 2 )even at high x. Introducing the following parameterization 

F 2N {x,Q 2 ) - A^F 2 (x,Q 2 ) (26) 

we obtain for the coefficient Pa results which presented in the Fig. [5] As it seems from the Fig[5]the coefficient 
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Figure 4: The a a coefficient of Ea W^ for different nuclei at different Q 2 . 
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Figure 5: The Rp 2 ration for different nuclei at different Q 2 . 
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Figure 6: The (3a coefficient of Ea [EE\ for different nuclei at different Q 2 . 
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Figure 7: Anomalous dimension 7 for the case of DIS on the proton at different Q 2 : Q 2 — 2.5 GeV 2 (solid 
line), Q 2 = 12 GeV 2 (dashed line), Q 2 = 60 GeV 2 (dotted line), Q 2 = 120 GeV 2 (dashed- dotted line). 

(3a is less sensitive to the details of the process, i.e. to the values of Q 2 and type of the nucleus than 
coefficient from FigJS] It is also interesting to note , that even at high values of x the coefficient (3 a does not 
equal to one, in contrary to the ola- 

5 Anomalous dimension 



Let's consider the definition of the average anomalous dimension 7 in DIS process via the integrated gluon 
density function 

xG(x,Q 2 ) oc (Q 2 ) 7 (27) 
see [36] . In this case the calculation of 7 is straightforward 

d\n(xG(x,Q 2 )) 



7 



(28) 



<91nQ 2 

The result for 7 for the case of DIS process on the proton is represented in FigfT] and for the case of DIS on 
the nuclei in the FigJSJ Considering the 7p anomalous dimension from the FigJT] it is interesting to note, 
that for DIS process on the proton at Q 2 > 60 GeV 2 the value of jp is below the value of BFKL anomalous 
dimension 'Jbfkl ~ 0.5 at whole range of x. At small values of Q 2 the 7p value is about BFKL 0.5 already 
at x oc 10~ 4 — 10 -5 , that one more time underline the importance of rescattering correction (shadowing 
corrections) in this kinematic region. From the result of the calculations of anomalous dimension for the 
different nuclei ^a in Fig(5] we see, that in contrary to the calculations of [23] the ja shows clear dependence 
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Figure 8: Anomalous dimension 7 for the case of DIS on the nuclei: A~20 (solid line), A=70 (dashed line), 
A=150 (dotted line), A — 197 (dashed- dotted line). 
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on the atomic number A at all values of Q 2 at small x < 10~ 4 . Decrease of the value of Q 2 in the DIS process 
leads to the increase of the value of 7,4, that indicates the increasing of value of the shadowing corrections in 
the process. 



6 Saturation momenta 

There are different definitions of the saturation momenta, which are used through the literature about the 
subject. For example, in papers |23) the saturation momenta in DIS process was defined with the help of a 
packing factor k p (Q 2 , x). As a saturation momenta there was considered a momenta where 

k p (Q%x) = 1/2 (29) 

at fixed x and impact parameter (if, of course, the impact parameter is introduced in definition of k p (Q 2 ,x)). 
In our paper we use a different definition of saturation momenta, borrowed from |37j . Following the definition 
of |37j we define a saturation momenta as a momenta where a maximum of the unintegrated gluon density 
function is reached at fixed impact parameter and fixed x: 

Q s (b,x) : f( X:< ^ S: ^ > f( x -> p & ) f or any k ■ k min < k < k max , (30) 

where k m i n and k max are correspondingly minimum and maximum momenta used in numerical calculations. 
This Q 2 s (b) definitely depends on impact parameter and, in fact, may be used in order to introduce the impact 
parameter dependence in the scaling solution of usual BK equation. In the next subsections we will use the 
definition Ea l30l for the calculations of saturation momenta of the proton and different nuclei. 



6.1 Saturation momenta for DIS process on the proton 

Plot of the saturation momenta of the proton, defined through the Eal3Ul are presented in the Fig.9-Fig[TU] . 
Considering the saturation momenta at fixed impact parameters, we could write a very simple expression for 
the approximate parameterization of the saturation momenta of the proton 

Q 2 (b,x) = Q 2 S0 (b) + Q 2 sl (b) ) d(b) , (31) 

where the coefficients Q| (&) , Q%i(b) and d(b) in this parameterization could be extracted from the Fig. 9 data. 
The fitting procedure for the data in the range of x such that x = 1.66 1CP 2 — 3.75 10 -8 gives the following 
values of coefficients 

Q 2 s (b,x) = F s (x)S(b) = f 8.69 + 12.65 —— J j GeV 2 , (32) 

with the coefficient d which does not dependent on impact parameter and with the proton radius from the 
Table [TJ Comparing the expression of the saturation momenta EqJ32]with the coefficient C from the expression 
Eq|9l we see, that the expression Eal32l gives C = 0.11 GeV 2 instead C — 0.0855 GeV 2 in initial conditions 
EqlHl The difference between these two values is due the "averaging" procedure used in fitting the data. The 
expression Eg J32I is the result of fitting of many points, and there is not necessary that the resulting curve will 
cross precisely the initial point in the fitted data. Obtained value of the coefficient d — 0.46 is close to the 
results of [37 . We obtained c = 2.2 instead c = 2.06 in the terms of the paper [57]. From the Eq. (|3"2"j) and 
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Figure 9: Saturation momenta of the proton as a junction of x and impact parameter. 




Figure 10: Saturation momenta as a junction oj impact parameter at fixed values of x in the left plot and the 
comparison of the obtained saturation momenta at b = with the GBW saturation momenta from the in 
the right plot of the picture. 
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Figure 11: Unintegrated gluon density function f(x, k 2 , b) / k 2 as a function of scaling variable r at b = and 
different x = 6.1 1CP 3 — 3.8 1CP 8 (curves from up to down in the left half of the plot) for the case of DIS on 
the proton. 

FigJTU]we see, that our answer for the saturation momenta is different from the obtained in the GBW model 
|17j . We will discuss these differences in the conclusion of the paper. 

With the help of Q|(6, x) it is easy to investigate scaling properties of the unintegrated gluon density function 
f(y,k 2 ,b). Plotting f(x,k 2 ,b) / k 2 as a function of only r = ^ at & = 0, we obtain results presented in 
the FigJTTJ As it seems from the FigfTI] the only approximate scaling of f(x, k 2 ,b) / k 2 exists, if we collect all 
the data at different values of x . This approximate scaling behavior is mostly pronounced near the maximum 
of the f(x, k 2 ,b) / k 2 function and it is clearly broken in the "tales" of the function at large and small values 
of r. Nevertheless, if we consider the f(x, k 2 ,b) / k 2 at fixed x and at different b, then we see a perfect scaling 
behavior of the function, see Fig U2l The only approximation scaling behavior of the unintegrated gluon density 
means, that f(x,k 2 ,b) / k 2 is likely a function of two variable, i.e. that f(x,k 2 ,b) / k 2 — F(x,r(x,b)) in the 
given framework. 



6.2 Saturation momenta for DIS process on the nuclei 

The behavior of the saturation momenta of the different nuclei as a function of impact parameter and x is 
presented in the FigfT3land in the left plot of the FigfT4l We parametrize the saturation momenta of the nuclei 
in the following form 

Q 2 SA (b,x) = Q 2 S0A (b) + Q 2 slA (b) (^) d(6) ■ (33) 
The fitting of the data gives the following expression for the saturation momenta of the nuclei 

Q 2 SA (b,x) = F SA {x)AS(b) = ^8.84 + 12.56 J AS ( b ) GeV 2 , (34) 

where S(b) is a Wood-Saxon nuclei profile function from EqQ31 We see, that the coefficients in the Fs(x) 
function in Eq 1341 are the same as in the expression for the saturation momenta of the proton in Eq l32l Indeed, 
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Figure 12: Unintegrated gluon density function f(x, k 2 ,b) / k 2 as a function of scaling variable r at fixed values 
of x and different b = — 12 GeV~ l for the case of DIS on the proton. 



IS 



Ne 



Zn 




Figure 13: Saturation momenta of nuclei as a junction of x and impact parameter. 



19 



if we compare the expression EqJMjwith the EqfT3]wc see, that the coefficient C from EqfT51must be the same as 
in EqJUJ Therefore, apart the impact parameter dependence and number of nucleons in target, the expressions 
for saturation momenta in Eqj32] and Eq[34] are the same. 

Now we could compare the expression Eq. (j3"2")l and Eq. It is easy to see the following relation between 

saturation scales 



ISA 



(0,x) 



Q 2 s (0,x)AS A (0) 
S p (0) 



(35) 



where S A (0) and S p (0) are the proton and nuclei profile functions correspondingly. So, from our expressions for 
the saturation momenta we obtain 

Q 2 SA (0,x)S p (0) 



Ql(0,x)AS A (0) 



1 



(36) 



The independent calculations of this ratio presented in the right plot of the Fig03] shows that the expression 
Eq. (1361) is indeed correct. Now it is easy to obtain a A n behavior of the nuclei saturation momenta. Having in 
mind that 



AS A (b) 
S p (b) 



= A l ' 3 k A {b) 



(37) 



we obtain 

Q 2 SA (b, x) = A 1 / 3 k A (b) Q%{b, x) . (38) 

We see, that we obtained a usual DGLAP A 1 / 3 dependence of the nuclei saturation momenta with the additional 
modification coefficient k A (b) which arises due the different sizes of proton and nuclei and which depend on the 
atomic number A of the nuclei. In the case of the considered nuclei this coefficient varies as k A (0) = 0.28 — 0.34. 
The DGLAP A 1 ' 3 dependence in Eq. <J3H]) 

again underline the relative smallness of the shadowing effects in the 
processes of interest. We also note, that k A (b) coefficient is related with the diluteness parameter k a from [2"T] 
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Figure 14: Normalized impact parameter profile of the saturation momenta of the proton and different nuclei in 
the left plot of the picture and ratio Ea \36\ for the different nuclei in the right plot of the picture. 

which is defined as 

n A {b) = A 1 ' 3 k A {b) . (39) 

For the considered nuclei this parameter varies as k a (Q) — 0.76 — 2 that support the point of view of [57] on 
the nucleus as on the pretty dilute system of nucleons. 



20 



Investigating the scaling properties of the nuclei unintcgratcd gluon density function in the same way as it 

.2 

was done before for the proton case and again introducing the r = — r variable we obtain results shown 
in the Fig. 15. As in the case of the DIS on the proton we see, that the scaling is precise only at fixed value of 




Figure 15: Unintegrated gluon density function f(x, k 2 , b) / k 2 as a function of scaling variable r at b = and 
different x — 6.1 10 -3 — 3.8 10 -8 (curves from up to down in the left half of the plot) for the case of DIS on 
the gold (A = 197). 

x. In general, the scaling behavior of the unintegrated gluon density function is only approximation, see again 
Fig. 15 and Fig. 16. 

7 Conclusion 

As the main result of this paper we consider the application of the BK evolution equation with the local impact 
parameter dependence to the DIS processes on the nuclei at small values of Bjorken x with the initial conditions 
of the rapidity evolution similar to the usual GBW ansatz. The precise form of the initial conditions for the 
case of DIS on the proton was obtained in the paper [4] with the help of F2 HERA data fit and we considered 
a DIS process on the nuclei using the same functional form of the initial conditions. Definitely, it is not fully 
clear why the main parameters of the initial conditions will not be changed if the processes with nuclei instead 
the proton are considered. Nevertheless, surprisingly, in the given framework we obtained the integrated gluon 
density function which similar to the integrated gluon density functions from the known parameterizations in 
the initial point of evolution, see Fig|2}Figi3J The explanation of this fact is very simple, in the given framework 
the conditions for the applicability of high energy formalism seems to be universal. Indeed, let's consider the a s 
coupling as a parameter which defines a physical "condition" of the process. We assume, that the value of a s 
is determined by some averaged saturation momenta. In considered range of x, x — 4.5 1CP 2 — 3.8 1CP 8 , the 
value of saturated momenta and, correspondingly the value of a s is almost the same for the proton and nuclei 
in LO calculation scheme. It gives a physical explanation of this universality, the application of the BFKL and 
BK equations in DIS process is independent on the considered target of the process. 

It was mentioned above, that the interesting property of the obtained results for the integrated gluon 
density function is that it stays in the range of the results for the integrated gluon density functions obtained 
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Figure 16: Unintegrated gluon density function f(x, k 2 ,b) / k 2 as a function of scaling variable r at fixed values 
of x and different b = — 110 GeV^ 1 for the case of DIS on the gold (A = 191). 
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iii the different parameterization schemes which base on the application of DGLAP equation to the low energy 
experimental data, see again FigJ3]-Fig|2] This fact, justifying chosen form of the initial conditions, causes 
a interesting question. In general, there is no physical reason why the curves obtained in framework of BK 
equation will be the same as the curves obtained in the DGLAP calculation scheme. Clearly, first of all, the 
reason for the such behavior of the curves is that the LO BFKL equation defines the same behavior for the 
integrated gluon density function as the LO DGLAP equation in some kinematic ranges. At large values of Q 2 , 
when BK triple Pomeron vertex corrections are small, the coincidence of the curves is a sequence of this fact. 
More interesting, that the curves obtained in BK equation formalism are still pretty close to the other curves 
at smaller values of Q 2 . As a explanation of this we assume that the triple Pomeron vertex corrections are 
relatively small in the considered kinematic range of x and Q 2 . As it seems, at given Q 2 and x the contribution 
of these corrections is not large, that makes resulting curves based on the BK and DGLAP equations formalism 
to be similar. Going into the region of smaller Q 2 at the small given and fixed values of x, we will see a larger 
deviations between the results. At very small values of Q 2 the considered LO BK equation fails to describe data, 
that could be a sign not only of a significance of the triple Pomeron vertex corrections but also a sign of the 
need to add a further, "net" diagrams corrections into the process, see [U [11] [T21 HS1 |38j [39] . Of course, going 
into the region of smaller x at fixed values of Q 2 , we also will obtain the more significant deviation between the 
curves of different approaches caused by the large triple Pomeron vertex corrections. 

The additional mark of the smallness of the triple Pomeron vertex corrections is a value of coefficient oia 
from the Eq[221 sec FigH] The LO DGLAP result gives a a — 1 for any value of Q 2 . As it seems from Fig [31 
even at Q 2 = 2.5 GeV 2 and x oc 10~ 3 the value of c*a ~ 0.8 for the gold target. The obtained value is close to 
one and it is larger then the value of a a obtained in [23] . for example. The same parameterization as in Eq|22l 
we can write for the F2N nuclei structure function, see Eq(55] and Eq|2l)l The calculations gives, that the F2N 
structure function is less sensitive to the shadowing, i.e. the coefficient Pa is closer to unity then coefficient cxa 
in the same kinematical region. 

The result of the anomalous dimension calculation for the case of the DIS on the proton, see Eg l28l and Figj7l 
shows a significant dependence of the 7p on the value of Q 2 . At large Q 2 and at small values of x, the 7p — > 1/2 , 
as it must be in the case of BFKL asymptotic results. For the smaller values of Q 2 the anomalous dimension 
is larger then half, that shows a corrections arise due the triple Pomeron vertex. The same calculations for 
the nuclei, see FiglHl show a dependence of the value anomalous dimension of the nuclei on the atomic number 
A. This dependence especially underlined for the case of the small values of Q 2 , in contrary to the results of 
calculations in |23j . 

Another problem, considered in the paper, is a problem of the saturation momenta in the DIS on the 
proton and on the nuclei. First of all, initially we assumed a similar and factorized form for the saturation 
momenta of the proton and of the nuclei in the process of interest, see Eq|9]and EafT3l The same factorized 
form was preserved in the final expressions for the saturation momenta of the proton and of the nuclei after 
the evolution over rapidity, see Eq[32] and EqIMl that in some sense justify the calculations of [531 [23] for 
example. Nevertheless, there is a principal difference between the calculations of the [35J [33] and present 
one which must be underlined. In spite of the introducing the form of impact parameter dependence in the 
solution of evolution equation after evolution, as it was done in 23, 25, 24 , the actual form of impact parameter 
dependence of the saturation momenta in the present framework is calculated through the rapidity evolution 
and, therefore the form of the saturation momenta with impact parameter dependence is different from the used 
in [231 [23 ED- 

The performed calculation of the saturation scale shows, that we did not observe a suppresion over atomic 
number in the expression for the saturation momenta, which was obtained in the similar calculations in different 



23 



papers, see for example [33] 123 EH S3 and references therein, see Eq. ([55]) and plot in FigfTJ] Still, we used a 
definition of saturation momenta different from the definition of [33] , but our results are also different from the 
results of [37] where the similar definition was used. The first simple fact, which could explain this result, is that 
we obtained and used a different expression for the saturation momenta. Our expression contains two parts, one 
part of this expression is a constant and second part is a function which grows with rapidity. This future of the 
considered model may be explained by the local impact parameter dependence introduced in the calculations. 
From the Fig.9-Fig ll3l we see, that in the large range of x the saturation momenta is almost a constant, the 
growth became at x oc 1CP 5 only. Considering the saturation momenta as the characteristic momenta of the 
scattering system which related with the averaged parton density of the system, we see, that until x oc 10 -5 
this density does not grow so fast. Therefore, the flatness of the saturation momenta as the function of x in 
broad small x region might be explained by the linear growth of target area in the impact parameter space 
accounted by the introduced impact parameter dependence. When a speed of increase of the area of the target 
is similar or larger than a speed of increase of the characteristic momenta then the saturation scale does not 
change so much. Only when the grow of the momenta is larger than the grow of the area of the target, only 
then we observe a increase of saturation momenta with the increasing of energy. We can conclude, therefore, 
that interplay between the linear growing of the area of the target in the impact parameter space and growing 
of the parton density "delays" growing of the saturation momenta. Only at asymptotically large energies, when 
the effect of the growth of the target area will be negligibly small, the expected suppression over the atomic 
number A perhaps could be observed . 

Relating the proton and nuclei saturation momenta, see expression Eq. (|38[) . we introduced a parameter kA(b) 
which is a parameter of the diluteness of the system, see 27 . As it was underlined in [37], see also references 
therein, this parameter is not so large for the case of real nucleus. Without discussing a new approaches 
introduced in [27] and a need of the additional rescattering corrections to the system amplitude, see for example 
[15] . we can conclude nevertheless, that the relative smallness of this parameter shows on the smallness of the 
triple vertex (shadowing) corrections as well. 

It is also interesting to compare obtained results with the other results on the saturation momenta. First of 
all, it is instructive to compare obtained form of the proton saturation momenta with the result of the GBW 
paper [T7], see FigfTUl As it seems from the plot, the compared saturation scales have different functional 
form and also, as a consequence, different numerical values at different values of x. At relative large values of 
x the obtained in this paper saturation scale is larger then the GBW one, whereas at very small values of x 
the GBW is larger. The explanation of this result was done before, the difference arises mainly due the very 
different functional form of the functions. The expression Eq. (|32p have a constant part which does depend on 
x, which arises because of the impact parameter dependence, and, therefore, the growth of the scale with the x 
decreases begins much later then in the GBW scale case. Therefore, it is difficult to compare our result with 
the different result based on the GBW scale parameterization and geometrical scaling application, see [44] for 
example. There the GBW saturation scale parameterization was used in order to extract the A n dependence of 
the nuclei saturation momenta. Obtained in these calculations coefficient n is different from DGLAP n = 1/3 
value obtained here, as it may be expected simply because of the different saturation scales parameterizations. 
Nevertheless, in the contrast to the results of [44] in the paper [43] the similar value of the coefficient n = 1/3 
for the nuclei saturation scale was obtained after the light nuclei data excluding from the fitting procedure. 
Still, it must be underlined, that the compared models have a different calculation frameworks, that makes they 
comparison pretty difficult, see Appendix A for some remarks on this subject. 

Obtained through the calculation expressions of the saturation momenta allow to investigate the scaling 
properties of the unintegrated gluon density function. In the given framework and in the considered kinematical 
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region we found that the scaling is only approximate future of the function as it seems from the Fig[TT]-Fig|T2] 
and Fig.15-Fig.16. In general the unintegrated gluon density function depends on the two variables, scaling 
variable r(b,x) and on the energy of the process (rapidity or value of x). Nevertheless, the scaling is precise 
when we fix the energy and consider different values of impact parameter. In this case, instead the different 
solutions at different impact parameters, the scaled solution arises and its depends only on the r(6) variable. 

Finally we would like to underline, that given calculations we consider as a first step in establishing of the 
framework for the investigation of more complex processes at very small values of x such as amplitude of proton- 
proton scattering at LNC energies or calculations of the amplitude of the exclusive Higgs boson production. 
The demonstration that in the framework of the interacting BFKL pomerons it is possible to describe a bulk 
of the DIS data on the proton and nuclei, leads us towards the application of the model in the calculations of 
these complex processes at the LHC collider energies. 
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Appendix A: 



In this appendix we would like to come back to some details concerning the DIS process on the proton from 
[3]. The main problem, which we faced in that calculations is a fail of the description of small Q 2 < 1 GeV 2 
DIS data for F% function, see details in 4, . If we will compare our approach with the saturation approach 
from [191 l29l [32] . which describe all data of DIS on the proton, we immediately will see that we missed a 
correct description of gluon density function at small values of momenta, i.e. we do not reproduce the initial 
condition of DGLAP evolution equation for the DIS process. This result may be explained by the absence of 
the additional "net" diagrams in the BK approach. Indeed, when in the DIS process Q 2 < Q 2 S then the use 
of only "fan" diagrams is acceptable. But when Q 2 > Q|, i.e. when the size of projectile is large then the 
averaged size of partons in the proton, then we definitely miss additional contributions in the amplitude which 
represented by the "net" diagrams, see |35]. Taking into account that we obtain for the proton saturation 
momenta Q 2 S « 1 GeV 2 we see that indeed we have to face a problems when Q 2 < 1 GeV 2 . 

Another hint to the solution of this problem is that during the calculations we obtain that the triple pomeron 
corrections are relatively small, see also Eq. ((39|) and [27]. Therefore, in spite of the long small x evolution we, as 
it seems, missed a correct initial DGLAP type condition in the region of small Q 2 S which is crucially important 
in this kinematic region. So, as a resolution of the problem, it will be interesting to combine the DGLAP type 
initial conditions from the |19j with the BK evolution equation with local impact parameter dependence. This 
task we will leave for another publications. 
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Figure 17: Initial profile of the unintegrated gluon density function f(x,k 2 ,b) as a function of momenta k 2 at 
fixed values of x and different values of b for the case of DIS on the proton.. 



As a illustration of the difference between the present model and model of [TH] we present here a plots of 
the unintegrated gluon density function Eq. ([3]) as a function of momenta at different values of x and impact 
parameter. This function has a dimension (as opposed to the dimensionless fix, k 2 ,b) / k 2 function), therefore 
the numerical values of the function in the plots is not important and we also skip the dimension of the function 
on the plots. Nevertheless, it is interesting to trace a creation of the large k 2 tale of this function beginning from 
the initial condition, FigJTl and through the small x evolution, FigjT] Comparing this plots with the similar ones 
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from the we see that our functions is absolutely different in the region of large k 2 , that may be explained 
as a result of the small x evolution. But there is a doubt that this tail explains a fail of the approach in the 
description of the F2 function in the region of small Q 2 . This large k 2 tail affects on the exclusive processes 
mainly, whereas the inclusive ones are sensitive to the region of small k 2 which seems to be not so different from 
the presented in the [T5] . 



x=0.136e-2 



0.05 



0.03 



0.02 



0.01 











b=0 GeV 
b=2.4 GeV 
b=4.8 GeV 


-1 * + * 
-1 

-1 „ 






































— v i y 




* 









0.001 



0.01 



0.1 



1 

k 2 [GeV 2 ] 



10 



100 



1000 



0.3 



0.15 



0.05 













b=0 GeV"] 
b=2.4 GeV"] 
b=4.8 GeV" * + 


































































x » 111 m -J 


«■» iici'l 




* * i * * 


* * * U * 









0.0001 0.001 0.01 0.1 1 10 100 1000 10000 

k 2 [GeV 2 ] 



x=0.124e-5 



x=0.375e-7 



3.5 
3 
2.5 

C 2 
^ 1.5 
1 

0.5 












b=0 GeV"] 
b=2.4 GeV, 
































































« iii iii iii r m in in « 


11 « in in hi IBM 









0.0001 



0.01 



1 

k 2 [GeV 2 ] 



100 



10000 



1e+06 



45 
40 
35 
30 
25 
20 
15 
10 
5 




0.0001 



0.01 











b=0GeV"] 
b=2.4 GeV" 












b=4.8 GeV 


1 X 






















































































***** 


iinonoonood 













100 
'■ [GeV 2 ] 



10000 



1e+06 



1e+08 
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